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ABSTRACT 

We present the results of a calculation of recombination coefficients for 0 2+ + e~ 
using an intermediate coupling treatment that fully accounts for the dependence of 
the distribution of population among the ground levels of 0 2+ on electron density and 
temperature. The calculation is extended down to low electron temperatures where 
dielectronic recombination arising from Rydberg states converging on the 0 2+ ground 
levels is an important process. The data, which consist of emission coefficients for 
8889 recombination lines and recombination coefficients for the ground and metastable 
states of 0 + are in Cases A, B and C, and are organised as a function of the electron 
temperature and number density, as well as wavelength. An interactive fortran 77 data 
server is also provided as an accessory for mining the line emission coefficients and 
obtaining Lagrange interpolated values for any choice of the two variables between 
the explicitly provided values for any set of wavelengths. Some illustrations of the 
application of the new data to nebular observations are also provided. 

Key words: atomic data - atomic processes - radiation mechanisms: general - 
radiation mechanisms: non-thermal - planetary nebulae: general - ISM: H II regions 
- astrochemistry - plasmas - methods: numerical. 


1 INTRODUCTION 

O II recombination lines have been observed in the spectra 


1993 

Liu et al 1995; Baldwin et al 2000 ; Garnett Diner- 

stein 

2001;|Wesson et a/|2003J Esteban et aS[2004| Wesson et 

alj2005 Robertson- Tessi & Garnett|2005; Peimbert & Pcim- 

bertJ2005||Liu et a/|2006; Esteban et al 

2009; 

Mesa-Delgado 

& Esteban]2010; Peimbert & Peimbert 

2013; 

Peimbert et al 


2014), and hence many of these lines (e.g. AA4076, 4089 and 


the VI multiplet) have been used ( 

Liu et al 

2004 

Tsamis 

et a|2004. Wesson et a|2005||Bastin:2006 Bastin & Storey 

2006| 

Wang & Liu 

2007 

McNabb et al 

2013| as diagnos- 


tic tools for probing the physical conditions of the emitting 
regions 

Regarding previous theoretical and computational work 


on 0 + recombination processes, Nussbaumer & Storey 


( 1984 1 calculated the effective dielectronic recombination co- 
efficients for a number of selected O II recombination lines 
in an LS-coupling approximation. |Pequignot et q/j ( |1991| ) 
calculated total and effective radiative recombination coef- 
ficients for several important O II lines and provided them 
in the form of fitting formulae valid over certain tempera- 
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ture ranges. Similar calculations were carried out by |Smits| 
( 1991 ) for selected lines where the calculations were based on 
an LS'-coupling scheme. The work of Nussbaumer fc Storeyl 


(1984) was extended by Storey (19941 who computed effec- 
tive recombination coefficients for a number of O II lines 
based on a full LS-coupled 0 + model atom at electron tem- 
peratures and number densities relevant to planetary nebu- 


lae. Liu et al (1995) provided O II radiative recombination 


coefficients for the 3d-3p and 4f-3d transitions based on the 


results of Storey ( 1994) but with transformation of radiative 
data to an intermediate coupling scheme limited to these 
transitions. All of the above calculations assumed that the 
three levels of the 0 2+ 3 P ground term are populated ac- 
cording to their statistical weights, which is not usually the 
case in the conditions typically found in photoionised nebu- 
lae. In general, for any ionic level, if the ambient density is 
below the critical density for that level, its population will 
depart from a Boltzmann distribution and the population of 
all states formed by recombination on to that level will be 
affected. The result is that for recombination on to an ion 
with more than one level in the ground term, such as O a+ , 
the relative intensities of lines within all recombination mul- 
tiplets will, in general, show some density dependence. 

R uiz et al 1 12003 ) investigated this effect empirically by 
considering the observed variation of the relative intensities 
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of the lines of 0 II multiplet VI as a function of the elec- 
tron number density derived from forbidden lines, mainly 


[Cl ill] . This approach was further refined by Peimbert et al 
( |2005[ ) and Peimbert & Peimbert (20051. We compare these 
authors’ results with our theoretical work in a subsequent 
section. Bastin ( 2006 1 and Bastin & Storey ( 2006 1 revised 


the calculation of Storey ( 1994 I by describing the whole re 


combination process in intermediate coupling and including 
the distribution of population among the 0 2+ levels. They 
used an R-matrix scattering code ( Berrington et al 1995 1 


to compute oscillator strengths and photoionisation cross- 
sections, and hence recombination coeffi cients, for the low 
lying states of 0 + . McNabb et al (20131 presented various 


O II diagnostic line ratios derived from the theoretical results 


of Bastin (2006) and used them to analyse a large number 
of planetary nebula and H II region spectra. 

In the present paper we use the techniques used by 


Bastin ( 2006 1 and 


Bastin & Storey (2006), and later also 


used by Fang et al (2011) for the recombination of N 2+ , but 
with some modifications and improvements as described be- 
low. The results presented here differ, in some cases signifi- 
cantly, from those obtained by Bastin (2006). We return to 


this point in the relevant subsequent sections. We provide 
the resulting data in a text file containing a list of transitions 
between the lower states of 0 + and emission coefficients for 
these lines as a function of electron temperature and den- 
sity and Case (A, B and C which will be explained later). 
We also provide a file containing recombination coefficients 
for the ground and metastable states of 0 + as functions 
of the same quantities. We also provide an interactive data 
server, in the form of a fortran 77 program, for mining the 
emission coefficient data and obtaining interpolated results 
for densities and temperatures between the explicitly com- 
puted values. More details will be given in the forthcoming 
sections. 


As discussed briefly in Storey & Sochi (2015b I, the cur- 


rent paper comes as part of a series of papers by the au- 


thors ( Storey fe Sochi||2013 2014 Storey et aZ|2014 Storey 


fc Sochi||2015a|b I aimed at providing the means to investi- 

gate the long-standing puzzle of the inconsistency between 
the results of elemental abundance and electron tempera- 
ture as derived from the optical recombination lines (ORL) 
with those derived from the collisionally-excited lines (CEL). 
Several possible explanations for the discrepancy have been 
proposed including the proposal that the recombination lines 
are formed in regions of much lower temperature than im- 
plied by the CEL diagnostics. We therefore extend the cal- 
culation of recombination coefficients down to much lower 
electron temperatures (100 K) than is usually the case. 

It has also been suggested that the free electron en- 
ergy distribution may not be best described by a Maxwell- 
Boltzmann distribution (Nicholls et al 2012 20131, so we 


plan to publish a further paper where recombination co- 
efficients have been calculated with the free electron ener- 
gies distributed according to the k. distribution suggested by 
theses authors. We note, however, that the non-Maxwellian 


model at present has little theoretical (Ferland et al 


2016) 


2016) 


or observational (Storey & Sochi 2014| Zhang et al 
support. 

The plan of the present paper is that in section[2]we give 
details of the calculation of the atomic parameters needed 
to calculate the population structure, including details of a 


new R-matrix calculation of radiative properties. In section[3] 
we give details about the structure and contents of the data 
files and provide general instructions and explanations about 
how the data should be explored and used. Section[4]presents 
some comparisons between theory and observation and the 
paper is concluded in section [5] with general conclusions. 


2 ATOMIC PARAMETERS FOR CT 


In this section we summarise the main features of the past 
and present calculations. We follow the methods used by 


Storey (19941, Bastin (2006), Bastin & Storey (2006) and 


Fang et al (2011 ) to compute the atomic parameters required 
to compute populations of the states of 0 + due to recombi- 
nation processes. We refer the reader to those references for 
a fuller explanation. 


2.1 The Term Scheme 


The principal series of 0 + is 2s 2 2p 2 ( 3 P)nZ which gives rise 
to doublet and quartet terms. Embedded within this series 
are the 2s 2 2p 2 ( 1 D)nZ, 2s 2 2p 2 ( 1 S)nZ, and 2s 2p 3 ( 5 S°)nZ series 
and terms of the 2s 2p 4 configuration. Some of the members 
of these additional configurations lie above the first ioni- 
sation limit and can be a source of population via dielec- 


tronic recombination. In the calculation of Storey (1994) 


the states and radiative transitions of 0 + were described 
in L5'-coupling but in this approximation several important 
process es cou ld not b e incorporated. To remedy this short- 
(2006) described the states of 0 + in a J c j; J 


Bastin 


coming, 

coupling scheme, where J c is the total angular momentum 
quantum number of the 0 2+ state, j is that of the valence 
electron and J is the total for the 0 + state. This enabled two 
significant extensions to the work of |Storey| ( 1994[ ) : firstly, 
the populations of the 3 P j c levels of the 0 2+ recombining ion 
could be explicitly accounted for, and secondly the three se- 
ries of Rydberg states converging on the 3 P j c levels could be 
treated separately. Figure [T] illustrates these states schemat- 
ically. These Rydberg states are a channel for a dielectronic 
recombination process that operates at very low tempera- 
tures, below approx imatel y 1000 _K. In their calculation of 
N 2+ recombination, Fang et al (2011) used the methods de- 
scribed by B astlnj ( |2006| ) but modified the coupling scheme 
describing the 0 + states to a J c l(K)s\ J scheme, where l 
and s = 1/2 are the orbital and spin angular momenta for 
the valence electron. In this scheme J c is coupled with l 
to give the quantum number K, which in turn is coupled 
with s to give the total angular momentum quantum num- 
ber J. This scheme has the advantage that it is the natural 
scheme for the high-Z Rydberg states which appear in en- 
ergetically close pairs characterised by J = K ± 1/2. It is 
also the scheme used to describe the continuum states in the 
R-matrix electron scattering codes described below in Sec- 
tion |2.3| Th e experimental resu lts for the energy levels of 
0 + given by Martin et al ( 1993 1 were used here to assign K 
and J quantum numbers to 0 + states with l > 2. For states 
of lower l a transition to LS coupling takes place which is 
fully accounted for in the intermediate coupling R-matrix 
treatment of the bound states described below. 
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Figure 1 . Schematic diagram of low-lying energy states of O 1 
with some low-temperature dielectronic recombination transitions 
between fine-structure autoionising levels below the three lowest 
ionisation thresholds of 0 2+ : ’ f’o , 3 Pi and 'T-fi. Some multiplets 
arising from transitions between low-lying terms are also depicted 
schematically. The figure is qualitatively informative but is not 
drawn to scale. 


2.2 0 + Populations 


Fang et al (2011) have described the calculation of the de- 


parture coefficients b(J c nl(K)s ; J) defined by: 


N(J c nl(K)s; J) _ / N(J c nl{K)s ; J) 
N e N+(S c L c -J c ) ~ \N e N+(S c L c -,J c ) 


b(J c nl(K)s ; J) (1) 


where N(J c nl(K)s; J), N + (S C L C -, J c ) and N e are the num- 
ber densities of the 0 + state J c nl(K)s; J, the 0 2+ state 
S C L C ; J c and free electrons respectively and the subscript 
S signifies the ratio of the number densities given by the 
Saha-Boltzmann equation. There are three stages to the cal- 
culation as described originally by Storey (] 1994| |. To sum- 
marise, the first stage involves a calculation of b(J c n) for all 
n ^ 1000 for the J c = 0, 1, 2 of the 3 P ground term of 0 2+ . 
In addition to the collisional-radiative processes described 
by |Hummer fc Storey| (1987 ) for hydrogen-like systems, au- 
toionisation probabilities to and dielectronic capture proba- 
bilities from the available continua, averaged over l for each 
n, are incorporated. For n > 1000 we assume b(J c n) = 1. 
In the second stage, departure coefficients b(J c nl(K)s ; J) 
are calculated, still in a hydrogen-like approximation for all 
n < m, where ni is the boundary above which (-changing 
collisions are sufficiently rapid that b(nl) = b(n ) is a good 
approximation. In the first two stages there are necessarily 
no radiative or collisional processes linking states of different 
J c - The third and final stage deals with the departure coeffi- 
cients for the states with n < rid, where rid is the boundary 
below which only radiative processes determine the popula- 
tions, i.e. where collisional processes can be neglected. 
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2.3 R-matrix Calculation for 0 2+ -f- e 


For the low-lying states of 0 + , the recombination coefficients 
were computed directly by integrating the photoionisation 
cross-sections, as obtained from an 0 2+ + e~ R-matrix cal- 
culation in an intermediate coupling scheme using the Breit- 
Pauli approximation. Cross-sections were obtained for all 
the states with principal quantum numbers n ^ 10 and or- 
bital angular momentum quantum numbers l ^ 4 in the 
principal series of 0 + plus all other states embedded among 
these. The calculation also provides bound state energies 
and oscillator strengths for use in the calculation of the 
0 + level populations. The main codes used for the calcu- 
lation are Autostructurj 2 (Eissner et q(||l974 Nussbaumer 


fe Storey! 1978| |Badnell|2011 ) which was used to create the 

atomic target, and the UCL-Belfast-Strath clyde version of 
the R-matrix codc0 ( Berrington et al 19951 to carry out the 
scattering calculations. 

The Autostructure code was used to generate the tar- 
get radial functions which are used as an input to the first 
stage of the R-matrix code. The radial data were produced 
using thirty-nine configurations synthesised from seven or- 
bitals; three spectroscopic (Is, 2s and 2p) and four correla- 
tion orbitals (3s, 3p, 3d and 4d). This configuration basis 
was previously used by |Storey et al (20141 in their work on 
the collision strengths of [O ill] forbidden lines and is given 
in Table (1) of ( Storey et a(| |2014| . The correlation orbitals 
are calculated in a Coulomb potential with central charge 
8 | \ni | • An iterative optimisation variational procedure was 
employed to obtain the orbital scaling parameters, A n i, as 


supplied in Table (2) of Storey et al (20141. 


In the R-matrix scattering calculations, 72 target terms, 


listed in Table (3) of Storey et al (20141, were used. For 
the 36 energetically lowest states, experimental energies ob- 


tained from Martin et al (1993), were used in the diagonal 


entries of the Hamiltonian in stage STG3 of the R-matrix 
code, instead of the theoretically computed values obtained 
from stage STG1, to provide accurate threshold energies and 
hence more accurate positions for resonances in the pho- 
toionisation cross-sections. 


2.4 Energy Levels 

The R-matrix calculation provides bound state energies for 
674 levels with principal quantum number n ^ 10 and 
orbital angular momentum quantum number l ^ 4. The 
Hamiltonian used in the calculation of bound state energies 
includes electrostatic terms and one-body relativistic terms, 
i.e. the mass and Darwin relativistic energy shifts and the 
spin-orbit interaction. Two-body spin and orbit dependent 
terms are not included in this formulation. In Table Q] we 
compare the calculated ionisation energies with experiment 
for all terms up to the highest of the 2s 2 2p 2 3d electron con- 
figuration, 26 terms in all. The Table also shows the total 
fine-structure splitting of each term compared to the exper- 
imental value. 


1 See Badnell: Autostructure write-up on WWW. URL: amdpp. 
phys . strath. ac .uk/autos/ver /WRITEUP 

See Badnell: R-matrix write-up on WWW. URL: amdpp. phys. 
strath . ac.uk/UK_RmaX/codes/ serial/WRITEUP 
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For states of the principal series, 2s 2 2p 2 ( 3 P)3 1 the agree- 
ment between theory and experiment is excellent with a 
maximum difference of 0.9% and differences that are gen- 
erally much less than 1%. There are larger differences for 
the members of the ( 1 D)nZ and series that are in 

this energy range, reaching a maximum of 1.1% and 2.8% 
respectively. For low l we can expect LS coupling to be 
a good approximation but to become increasingly inappro- 
priate as l increases, independent of n. Liu et al (1995) dis- 
cuss the coupling schemes that are appropriate for 0 + states 
and conclude that there is a partial breakdown of LS cou- 
pling for some of the ( 3 P)3d states and that LS coupling is 
completely inappropriate for the ( 3 P)4f levels. Departures 
from LS-coupling within a configuration ( 3 P )nl depend, to 
a first approximation, on the inverse of the energy separa- 
tion between states of the same total angular momentum 
J belonging to different terms, and the magnitude of the 
spin-orbit parameter for the orbital nl. The accuracy of the 
calculated spin-orbit parameter can be estimated by com- 
paring the total fine-structure splitting of individual terms 
with experiment. Table [l] shows that the calculated and ex- 
perimental total fine-structure splittings of the ( 3 P)3Z terms 
agree within 10%. 

The large differences between theory and experiment 
for the very small fine-structure splittings of the 2s 2 2p 3 2 D° 
and 2 P° terms arises because the spin-orbit parameter is 
zero for the half-filled 2p shell and the observed level sepa- 
rations are due to two-body fine-structure terms which are 
not included in our R-matrix formulation. In Table [2] we 
compare calculated and experimental energies for levels of 
the 2s 2 2p 2 ( 3 P)4f and 5f configurations. These levels are no 
longer well represented in LS-coupling but appear as groups 
attached to each of the 3 Pj parent levels. We show the sta- 
tistically weighted mean energies of each group and the total 
fine-structure splitting of the group. 


2.5 Bound-Bound and Bound-Free Radiative Data 


The R-matrix calculation provides electric dipole oscillator 
strengths between all 674 bound states, which include spin- 
changing intercombination transitions. It also provides pho- 
toionisation cross-sections resolved by final target level from 
which recombination coefficients can be calculated. The en- 
ergy mesh used to calculate the cross-sections is divided 
into two sections. From the energy of the 2s 2 2p 2 J P 2 tar- 
get state to just below the 2s 2 2p 2 4 D 2 threshold, a mesh 
of variable step is used. The step length is varied to ensure 
that resonance features in the cross-section, whose positions 
and widths have been previously determined, are well re- 
solved. Details of this approach can be found in |Kisielius ef| 
[a/ (1998). This mesh terminates at effective quantum num- 
ber, v = 10 relative to the 4 D 2 threshold and comprises 
14244 energies. For higher energies the energy mesh between 
thresholds is defined by a fixed step, A v = 0.001 in effec- 
tive quantum number relative to the next threshold, giving 
rise to a further 11569 energies. For v ^ 10 relative to any 
threshold, the cross-section is obtained by Gailitis averag- 
ing with the mesh determined by Av = 0.001 relative to 
the next highest threshold. The mesh extends to just above 
the 2s 2p 3 3 D° target state to an energy of 1.095 Rydberg. 
The calculations described by |Bastin| ( |2006| ) and |Bastin fcj 
Storey ( 2006 1 were also based on an R-matrix calculation of 


energy levels, oscillator strengths and photoionisation cross- 
sections, albeit with a simpler scattering target than the one 
used here. However, a re-examination of those photoionisa- 
tion cross-sections shows that some of the cross-sections for 
the more highly excited states display non-physical features, 
such as unexpected oscillations and elevated continua. The 
results of the present calculation are therefore to be pre- 
ferred. 

For the 4f-3d and 3d-3p transition arrays we can com- 
pare our energies and branching ratios with those of |Liu| 


et al ( 1995 1 who used the atomic structure code SUPER- 


STRUCTURE (Eissner et al 1974 Nussbaumer & Storey 


1978 1 to calculate energy levels and dipole matrix elements 


for these transition arrays in intermediate coupling. |Liu et al\ 
(1995) made empirical corrections to the term and level en- 
ergies and spin-orbit parameters to bring the energies close 
to experiment. In the R-matrix approach used here we can 
make corrections to the target energies but not to the final 
0 + bound state energies. In Table[3]we compare individual 
level energies with experiment and with the results of |Liu et] 


al ( 1995 I for the 3d and 4f configurations. The empirical ad- 


justments made by Liu et al ( 1995) mean that they obtained 
better fine-structure splittings than those obtained here. In 
Tables [4] and [5] we compare the branching ratios for selected 
stronger components of the 3d-3p and 4f-3d transition ar- 
rays computed by Liu et al ( 1995 I with the present work. 
For the larger branching ratios, greater than 0.1, there is 
general agreement between our results and those of |Liu et al\ 
( 1995 1 to within 3% for the 3d-3p transitions and to within 
10% for the 4f-3d transitions. However, Liu et al ( 1995 1 point 
out that whereas the 4 F, 2 P and 2 D terms of the 2s"2p 2 3d 
configuration are well described by L5'-coupling, the 4 P, 4 D 
and 2 F terms are not. In Table pi we see that the largest 
differe nces between the present work and that of |Liu et al\ 
( 1995 1 occurs for transitions involving the 4 D and 2 F terms 
and in particular the 4 D 5 / 2 and 2 F 5 / 2 levels where the dif- 
ferences reach 40%. It is probable that the branching ratios 
calculated by Liu et al ( 1995 I with empirical corrections to 
level energies and spin-orbit parameters are closer to real- 
ity than our ab initio results for these particular transitions. 
We note, however, that the predicted strength of the corre- 
sponding recombination lines also depends upon the relative 
populations of the 0 2+ 3 Pj ground levels and these were as- 
suinedto be in proportion to statistical weights by Li u et\ 
\gj ( jl995[ ) . In the present calculation the 0 2+ 3 Pj level pop- 
ulations were explicitly included and accounted for in the 
calculation of all the 0 + level populations. We return to 
this point below. 


2.6 Autoionisation Probabilities 

The calculation of the departure coefficients above the 3 Po 
ionisation limit and below the 3 P 2 limit requires, in addi- 
tion to the usual collisional-radiative processes, rates for au- 
toionisation and dielectronic capture. As shown in Figure [l] 
these additional processes apply for n ^ 63 in the ( 3 Pi )nl 
series and m ^ 38 in the ( 3 P 2 )ml series. However, in an 
R-matrix calculation there is no allowance for the competi- 
tion between autoionisation and radiative decay of the va- 
lence electron, so a simple integration of the photoionisation 
cross-section does not correctly treat the dielectronic recom- 
bination via these high n, m states. We therefore use AU- 
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TOSTRUCTURE (Eissner et al 1974 Nussbaumer & Storey 


1978 Badnell 2011) to calculate the autoionisation proba- 


bilities for selected n,m < 1000 and l ^ 30 and to avoid du- 
plication of the dielectronic component due to the ( 3 Pi )nl 
and ( 3 P 2 )ml Rydberg states we begin the calculation of the 
R-matrix photoionisation cross-sections from the low-lying 
states at the energy of the 0 2+ 3 P 2 level. This provides the 
dielectronic contribution to the recombination, allowing for 
all the usual collisional-radiative processes but not the non- 
resonant or radiative recombination component, which we 
obtain by extrapolation of only the background photoion- 
isation cross-section from above the 3 P 2 threshold to the 
energy region between the 3 Po and 3 P 2 states. 


3 RESULTS AND DISCUSSION 

3.1 Structure and Contents of Data File 

The O II recombination lines and their theoretically- 
computed emission coefficients (erg cm 3 s _1 per 0 2+ ion 
per electron) up to n = 7 are archived in a single file called 
“OIIlines_ABC” , where ABC refers to the three “Cases” for 
which emission coefficients were calculated. Case A refers 
to the situation where all 0 + radiative transitions are as- 
sumed optically thin. In Case B we assume that all radiative 
transitions terminating on the 0 + 4 S° /2 ground level are 
sufficiently optically thick that the photons are re-absorbed 
on-the-spot in the same transition. Hence, these decays are 
excluded from the calculation of level populations. For Case 
C we also exclude decays to the 0 + 2 D° /2 and 2 D° /2 . This 
is a simple and relatively crude treatment of the effects of 
optical depth in resonance transitions but comparison with 
observations indicates that Case B is a much better approx- 
imation to reality in the majority of PNe than Case A or 
Case C. We briefly discuss the possible relevance of Case C 
in Section |U 

The “OIIlines_ABC” Hie consists of two main sections 
preceded by a short header which consists of a general de- 
scription of the file structure and the abbreviations used to 
label and explain the data blocks and their units. The struc- 
ture of these sections is as follows: 

1. The first section consists of 8889 indexed text lines corre- 
sponding to the provided 8889 O II transitions where each 
text line contains full identification of the corresponding 
transition by its upper and lower level as well as its wave- 
length in A and other relevant atomic designations. 

2. The second section consists of 8889 parts where each part 
corresponds to one transition as indexed above. Each part 
is made of three blocks where each block corresponds to one 
of the three Cases. Each one of these blocks is a rectangular 
array of emission coefficients of the indexed transition ar- 
ranged in 25 rows, as a function of log 10 (T e [I<]) between 
2. 0-4. 4 in steps of 0.1, and 16 columns as a function of 
log 10 (A e [cm -3 ]) between 2. 0-5.0 in steps of 0.2. 

More details about the structure and contents of the data file 
are provided in the ReadMe file associating the distributed 
data set. 

The data set distributed as part of the current work also 
includes a file called “OIImeta^VBC” containing the recom- 
bination coefficients (cm 3 s -1 per O 2 " 1 " ion per electron) to 
the ground and meta-stable levels of 0 + as functions of elec- 
tron temperature and number density and Case (A, B and 


C). The file contains 14 data blocks. The structure of the 
file is as follows: 

1. A block for the effective recombination coefficients of the 
ground level 2s 2 2p 3 4 S° /2 in Case A. 

2. Four blocks for the effective recombination coefficients of 
the two levels 2s 2 2p 3 2 D° and 2 D° each in Case A and 

B. 

3. Six blocks for the effective recombination coefficients of 
the two levels 2s 2 2p 3 2 P° /2 and 2 P° /2 each in Case A, B and 

C. 

4. Three blocks for the total recombination coefficients, one 
block for each Case. 

Each data block is structured as a function of electron tem- 
perature and number density in log 10 identical to the struc- 
ture of the emission coefficient blocks in “OIIlines_ABC” 
file, as described above. We note that direct recombination 
and all cascade contributions, excluding those between the 
metastable states themselves, are included in the calculation 
of the coefficients in the “OIImeta_ABC” file. 


3.2 Data Mining 

The data in the “OIIlines_ABC” file can be excavated using 
the supplied fortran 77 interactive data server. This server 
reads the file and provides easy access to the lines and their 
emission coefficients at the desired T e and N e . When the 
required data are within the T e and N e ranges but are not 
at the grid points of T e and N e as described above, a 6-point 
Lagrange interpolation in temperature and density is used 
to provide the desired data. Two main options are offered 
by the data server: (1) Generate a line list, and (2) Extract 
emission coefficients for a user supplied list of lines in user 
specified conditions. 

As for the first option, the server offers two choices: (la) 
List all the 8889 lines (between 367.97 and 3.90713xl0 7 A), 
and (lb) Specify a wavelength range for the list of lines de- 
sired by the user. In each one of these cases (i.e. la and lb) 
two files are generated that contain the list of lines where the 
lines in one of these files (named “OIIdataJist_eorder”) are 
ordered in decreasing emission coefficient and in the other 
file (named “OIIdata_list_worder”) in increasing order of 
wavelength; the full list in the latter file is preceded by a list 
of only the strongest lines in the required wavelength range. 
The emission coefficients in these files correspond to typical 
nebular conditions of T e = 10 4 K and N e = 10 4 cm -3 . 

Regarding the second option, the user is asked to sup- 
ply a list of lines in a text file for which the data are re- 
quired. The user then has the choice between (2a) Providing 
a list of T e [K]-A e [cm -3 ] pairs in a file for which the data 
are needed, or (2b) Providing a range of log 10 (T e [K] ) and 
log 10 (A r e [cm -3 ]) in a file defining a one- or two-dimensional 
grid for which the data are required. In all of these options 
(i.e. la, lb, 2a and 2b) the user has the choice to have the 
data for Case A, B or C. More details about the data server 
and its input and output files are provided in the ReadMe 
file associating the distributed data and code. 
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Figure 2. Computed fractional intensities of the eight transitions 
of the VI multiplet of O II at T e = 10000 K and in Case B as a 
function of log| ( j[AC(crrr 3 )]. 


Figure 3. Computed fractional intensities of the eight transitions 
of the VI multiplet of the O II at T e = 1000 K and in Case B as 
a function of logjgfiVefcm -3 )]. 


4 EXAMPLES OF COMPARISONS WITH 
OBSERVATIONS 

4.1 Multiplet VI 

The lines of O II multiplet VI are among the brightest and 
most accessible spectroscopically. The relative intensities of 
the components of VI depend upon the electron density via 
the distribution of population in the 3 Pj ground levels of 
0 2+ . This is illustrated in Figure [ 2 ] which shows the frac- 
tional intensities of the eight lines of VI at a temperature of 
10000 K. The relative populations of the 3 Pj levels are ap- 
proaching their local thermodynamic equilibrium values at 
N e = 10 5 cm -3 . At the higher densities the strongest com- 
ponent is the 4649.13 A line which, arising from the 4 D° /2 
level, is primarily formed by recombination from the 0 2+ 
3 P 2 level at this temperature. The relative intensity of this 
component falls rapidly as the density an d the po pul ation 
of the 3 P 2 level fall. In their Figure 1, McNabb et al (20131 


1 2004 ) for the higher ionisation species such as 0 2+ in M42 


also present the fractional intensities of the eight lines of VI 
over the same range of densities. Our results, which super- 
sede those of McNabb et al (20131, are nevertheless in good 
agreement for this case. 

In Figure [3] we show the fractional intensities for multi- 
plet VI at 1000 K. At this temperature the density variation 
of 4649.13 A is much less marked than at 10000 K partly 
due to the effect of dielectronic recombination through the 
high-n Rydberg states converging on the 3 Pi and 3 P 2 lev- 
els. Through this mechanism ( 3 P 2 )rd levels can be populated 
even at low electron density and, by radiative decays, even- 
tually populate the upper level of the 4649.13 A line. This 
process is only significant when the mean thermal energy of 
the free electrons is comparable to the 3 Pj level separations. 


Very similar results were shown by McNabb et al (20131 in 
their Figure 1. 


Baldwin et al ( 2000 1 and Esteban et al ( 2004 1 have 
published deep optical spectra of the Orion Nebula, M42. 
If we assume an electron temperature we can use the mea- 
sured and computed relative intensities of the VI lines to 
determine an electron density for the emitting regions. We 
adopt the temperature of 8320 K derived by|Esteban et al\ 


and locate the electron density that gives the best fit (least 
squares) between theory and observation for the seven lines 
of multiplet VI for which measured intensities are available. 
In Figure[4]we show the theoretical fractional intensities plus 


the observed values from Baldwin et al ( 2000 1 and Esteban 
| et aZ| ( [2004] ) plotted at the best fit electron densities. The ver- 
tical error bars are the error estimates given by the authors 
and the horizontal error bar attached to the 4649.13 A line is 
the one sigma uncertainty in the density derived from those 
observational errors. The agreement between theory and ob- 
servation at the best fit densities is excellent. The derived 
densities of 7330 cm -3 and 4400 cm -3 can be com p ared 
with the value of 8900 cm -3 adopted for M42 by Esteban et 


adopted for M42 by 

[a/ ( 2004 1 from a weighted average of forbidden line density 
diagnostics. 

In Figure [5] we again use the observations of Baldwin 
et al ( 2000|) and |Esteban et al (2004) to derive electron 


densities which are the best fit to the seven observed lines 
of multiplet V2. The lines of V2 are weaker than those of 
VI and the observational errors are consequently relatively 
larger, leading to larger uncertainties in the derived den- 
sities. Within the error bars there is reasonable agreement 
between the results derived from multiplets VI and V2. 

In Figure[6]we plot the intensity of VI A4649.13 relative 
to the total VI intensity as a function of electron number 
density. The results from the current work (CW) were cal- 
culated in Case B at 10000 K. We also show the effect on 
the ratio of switching off fine-structure dielectronic recombi- 
nation (CW-NoFSDR). The solid curve (PP) shows the em- 
pirical results derived by Peimbert & Peimbert ( 2005 I from 


the relationship between the observed values for a range of 
objects and the density derived from forbidden line ratios, 
in practice mainly from [Cl ill]. Figure [6] also shows the re- 
sults we obtain from the best fit of the relative intensities of 
all the observed lines of multiplet VI to the theoretical rel- 
ative intensities for three types of object: H 11 regions, PNe 
with low ADF, which we define as ADF ^ 2.5, and PNe 
with high ADF > 2.5, w here the ADF values were taken 
from Tsamis et al (20041, Table 9 (ORL/OPT for 0 2+ ). 
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Figure 4. Computed fractional intensities in Case B of seven 
transitions of the VI multiplet of the O II recombination lines at 
T e = 8320 K as a function of log 10 [./Ve(cm -3 )] with the observa- 
tional intensities of the H II region Orion Nebula (M42) plotted 


at N e = 7330l|gQ cur 


(red right) as obtained from 

i—3 


Baldwin 


et a/||2000) and at N e = 4400^^40 cm 3 (blue left) as obtained 
ijE; 


trom ]Esteban et al\ ( |2004j ) . Some points are displaced horizontally 
to avoid overlapping. 



Figure 5. Computed fractional intensities in Case B of seven 


transitions of the V2 multiplet of the O II recombination lines at 
T e = 8320 K as a function of log 10 [Ar e (cm -3 )] with the observa- 


tional intensities of the H II region Orion Nebula (M42) plotted 
at N e = 10210^4g®Q 0 cm 


Baldwin 


3 (red right) as obtained from 

et al (2000) and at N e = 2170^,!,|q° cm -3 (blue left) as obtained 


trom jEsteban et q/| ( |2004[ . Some points are displaced horizontally 
to avoid overlapping. 


A4649.13 relative intensity with density is less pronounced 
at 1000 K than at 10000 K so that a given observed ratio 
would generally imply a lower electron density than those 
obtained from the fits at 10000 K and listed in Table [6] 
We note that the empirical curve of Peim bert fe Peimbert| 
1 2005 1 is in good agreement with our theory except for low 


density objects, where the [Cl ill] densities used by Peim- 


bert & Peimbert ( 2005 1 tend to be lower than those derived 


here from the recombination line theory. There is therefore 
no clear indication from Figure [6] that the recombination 
lines are formed in a region of significantly different density 
to the [Cl ill] forbidden lines, as might be expected if they 
originated from different physical regions. 

The density variations of the VI lines and the lines of 
other O II recombination multiplets result primarily from 
the density variations of the populations of the 0 2+ 3 Pj 
levels. The critical density Ni for a given level i, is given 
by: 


fc<i k^i 


(2) 


where A;*, are radiative transition probabilities and qik are 
rate coefficients for collisional excitation or de-excitation. If 
the ambient density is much greater than Ni for a given 
level, the level population will approach the value given by 
the Boltzmann distribution. At 10 4 K the critical densi- 
ties of the 3 Pi and 3 P 2 levels are approximately 610 and 
4350 cm -3 respectively but for the level populations to ap- 
proach within 10% of Boltzmann values will require densi- 
ties approximately 10 times the critical densities, i.e. 6100 
and 43500 cm -3 . Typical PN (and H II region) densities 
are lower than this so we can conclude that the lines of all 
O II recombination multiplets will normally exhibit a signif- 
icant dependence on density. Thus only in very high density 
nebulae can it be assumed that the 0 2+ 3 P,j levels have a 
Boltzmann distribution of population and that the previ- 
ous LS-coupling theoretical results (e.g. |Storey| ( jl994[ >) can 
be used to infer abundances from individual lines without 
knowledge of the electron density. 

An additional consideration is that critical densities de- 
crease with decreasing temperature, so that at 10 3 K for ex- 
ample, the values of 10 x AT rlt for the 3 Pi and 3 P 2 levels of 
0 2+ are 2500 and 15500 cm -3 respectively. These values are 
significantly lower than those at 10 4 K but not sufficiently 
low that a Boltzmann distribution can be safely assumed. 

It is useful to apply the same argument to the recombi- 
nation lines of N II and Ne II, where we find that 10 x Nf rlt 
is 20500 cm -3 for the 2 P° /2 level of N 2+ and 1.96 x 10° 
and 2.81 x 10 s cm -3 for the 3 Pi and 3 P 0 levels of Ne 2+ . 
Therefore we anticipate that all recombination multiplets of 
N II and Ne II will show significant density dependence at 
nebular densities. 


The densities resulting from the fits are given in numerical 
form in Table [6] which also lists the [O ill] optical forbid- 
den line temperature for each object. The theoretical values 
used in the fitting procedure were all calculated at 10000 K, 
which may not be appropriate for the high ADF PNe where 
there is some evidence that the recombination lines origi- 
nate from a region at significantly lower temperature. As 
discussed above and shown in Figure [3] the variation of the 


4.2 Density and Temperature Determinations 

In the previous section we chose to use observations of 
M42 to initially illustrate the determination of electron 
density because there is broad agreement about the elec- 
tron temperatures in M42 and H II regions generally. The 
same cannot be said about PNe where some models pro- 
pose that recombination lines are formed in regions of much 
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Figure 6. The intensity of A4649.13 relative to the total multiplet 
VI intensity as a function of log 1 Q[Al e (cm — 3 )] from [Peimbert fc| 
|Peimbert| ( |2005[ ) (PP), the current work (CW) and the current 
work neglecting fine-structure dielectronic recombination (CW- 
NoFSDR). The two curves showing the present work correspond 
to T e = 10000 K in Case B. See Table [6] for information about 
the observational points. 


lower temperature than those derived from forbidden line 
diagnostics. It is therefore of interest to seek diagnostics 
that determine both the density and the temperature of 
the regions where recombination lines are formed. In Fig- 
ure [7] we plot a predominantly density sensitive line ratio 
A4649.13/A4661.63 against one which is mainly sensitive to 
temperature A4649.13/A4089.29, in Case B. The former pair 
of lines are from the VI multiplet and the latter pair com- 
prise the ratio of a VI line with the strongest line from 
the 4f-3d transition array, A4089.29. In the latter case, both 
lines arise from the highest J of the transition array and 
hence depend strongly on the population of the 3 P 2 ground 
level of 0 2+ , so that there is relatively little sensitivity to 
density in the ratio. Figure 2 of McNabb et al (2013) plots 
the same pair of ratios based on the earlier calculations of 
Bastin (2006| ), which shows some minor differences to our 
Figure [7J particularly in the values of the temperature sen- 
sitive A4649.13/A4089.29 ratio. As stated above, the earlier 
results contained defective photoionisation cross-sections for 
some states and are therefore superseded by the present 
work. 


Peimbert & Peimbert (20131 have suggested that, in low 


resolution spectra, the intensity of the A4089.29 line may be 
affected by blending with a Si IV line at 4088. 86A which, 
once corrected for, would potentially increase the observed 
A4649.13/A4089.29 ratio and result in deriving higher elec- 
tron temperatures from this ratio. The Si IV A4088.86 line 
is part of a doublet with the other, weaker component at 
A4116.10, and a line at this wavelength has indeed been 


observed in M42 and several PNe (Peimbert & Peimbert 


2013), although other identifications have been proposed. 


Thus caution should be exercised in deriving electron tem- 
peratures using a single ratio such as A4649.13/A4089.29. 

Notwithstanding the possible blending issues, we show 
observational values for 23 objects including H II regions, 
low ADF PNe and high ADF PNe, defined as in the previous 


section. The sources of the observational data are as follows. 


H II regions: Baldwin et al (20001, Tsamis et al (2003), Este- 

ban et al (2004), Garcia-Rojas et al 

12004). Low ADF PNe: 

Hyung et al (2000), Tsamis et al 

2003), Garcia- Rojas et al 

(20151. High ADF PNe: Liu et al 

>. 

19951, Liu et al 

(2000), 

Liu et al ( 2001j) , Tsamis et at (2003 

Liu et al (2006), 

Fang & 

Liu (2011), McNabb et al (2016). The A4089.29 line which is 


used to constrain the temperature is a 4f-3d transition and 
is relatively weak. This is reflected in the large error bars 
on the A4649.13/A4089.29 ratio particularly for the H II re- 
gions while in the high ADF PNe, which exhibit stronger 
recombination lines, the error bars are relatively smaller. 
The three classes of objects broadly occupy different regions 
of the diagram. Most marked is the tendency for the high 
ADF PNe to cluster in the very low temperature region of 
the plot. The low ADF PNe, on the other hand, mainly oc- 
cupy the region where temperatures are closer to 10000 K, 
although there are exceptions. The data for H II regions show 
a surprisingly large scatter in terms of temperature, which 
may be partially attributable to the intrinsic weakness of 
recombination lines when ADF factors are close to unity. 
We emphasise also that the theoretical results in Figure [7] 
were calculated in Case B which may not be the best ap- 
proximation for H II regions which have large spatial extent 
and low ionisation. This may lead to sufficient population to 
build up in the 2 D° and 2 D° levels of 0 + for radiative 
transitions to these two levels to become optically thick and 
hence for Case C to be more appropriate than Case B. 


5 CONCLUSIONS 

We have calculated emission coefficients for recombination 
lines of O II over a range of electron temperatures and den- 
sities in Cases A, B and C, suitable for application to the 
spectra of photoionised nebulae. The calculation treats the 
important states in intermediate coupling and uses rates for 
bound-bound and bound-free radiative processes computed 
with the R-matrix method of electron scattering. The pop- 
ulations of the fine-structure levels of the ground term of 
0 2+ are explicitly included, as is the effect of dielectronic re- 
combination via high-n Rydberg states converging on these 
levels. To the best of our knowledge, these processes have 
not been included in any previous calculation of recombina- 
tion coefficients for O II. Based on the above calculations, a 
large data set that comprises 8889 recombination fines with 
their emission coefficients as functions of electron tempera- 
ture, number density and Case (A, B and C) accompanies 
this paper. An interactive data server, in the form of a for- 
tran 77 code, is provided for the convenience of the users 
to explore the data and obtain Lagrange interpolated val- 
ues between the explicitly provided ones. We also provide 
a data set of coefficients for recombination to the ground 
level and the four metastable levels of 0 + plus the total 
0 2+ + e“ recombination coefficients as functions of electron 
temperature and density and Case, for use in modelling the 
oxygen ionisation balance and calculating the recombination 
contribution to the excitation of the O II forbidden lines. 

Comparison with observations showed that the lines of 
multiplet VI, the brightest of the O II recombination spec- 
trum, provide a novel means of determining electron density, 
at least in those objects where the electron temperature is 
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Figure 7. The ratio of our theoretical line intensities 7(A4649.13)//(A4661.63) versus the ratio of our theoretical line intensities 
7(A4649.13)//(A4089.29) in Case B and as functions of log 10 [T e (K)] (solid red) and log 10 [A r e (cm -3 )] (dashed blue). See the text for 
the sources of the observational points. 


reasonably well known. We have also shown that there is 
good agreement between our calculated relative strengths of 
multiplet VI as a function of electron density and empirical 
results derived from observations and CEL densities. The 
density dependence of the lines of the V 1 multiplet at typi- 
cal nebular densities is a result of these densities being below 
the critical densities for the ground levels of 0 2+ , and sim- 
ilar density dependence is therefore expected to appear in 
the components of all recombination multiplets of O II and 
in all ions where a similar situation prevails. Determining 
electron temperature from O II recombination lines is more 
problematic due to the relative weakness of the 4f-3d lines 
which, when compared to multiplet VI, show sensitivity to 
temperature. There does, however, seem to be evidence that, 
for the group of high ADF PNe considered here, the recom- 
bination lines are formed at a markedly lower temperature 
than is derived from CELs. 
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7 STATEMENT 

The complete data generated in this work can be obtained 
in electronic format with full precision from the Centre de 
Donnees astronomiques de Strasbourg (CDS) database. The 
fortran 77 data server is compiled and tested thoroughly 
using gfortran, f77, and Intel fortran compilers on Ubuntu 
12.04 and Scientific Linux platforms. Representative sample 
results from all these compilers and on all those platforms 
are compared and found to be identical within the stated 
numerical accuracy. 
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Table 1. The 26 lowest bound terms of the 0 + ion, which include all levels up to the highest of the ls 2 2s 2 2p 2 ( 3 P)3d configuration in 
experimental energy order, and their experimental ( E eyi ) and theoretical (E t h) energies in Rydberg relative to the first ionisation limit 
(ls 2 2s 2 2p 2 3 P Q ), as well as the total fine-structure splitting in cm -1 of each term experimentally (TSTE) and theoretically (TSTT). The 

column %E represents the percentage relative difference between the experimental and theoretical energies 10Q (^ X — ^ihlj w hii e the %S 
column represents a similar percentage for the splitting. The experimental energies are obtained from the NIST database (www.nist .gov) 
while the theoretical energies are obtained from the bound stage of the R-matrix code. The Is 2 core is suppressed from all configurations. 


Index 

Term 

f'ex 

Eth 

%E 

TSTE 

TSTT 

%S 

i 

2s 2 2p 3 4 S° 

-2.581409 

-2.587200 

-0.22 

0.0 

0.0 

0.0 

2 

2s 2 2p 3 2 D° 

-2.337020 

-2.333770 

0.14 

20.0 

00 

bo 

143.9 

3 

2s 2 2p 3 2 P° 

-2.212625 

-2.203115 

0.43 

2.0 

-6.6 

430.9 

4 

2s 2p 4 4 P 

-1.488604 

-1.487096 

0.10 

245.6 

228.2 

7.1 

5 

2s 2p 4 2 D 

-1.068751 

-1.052943 

1.48 

8.0 

-6.6 

181.9 

6 

2s 2 2p 2 ( 3 P)3s 4 P 

-0.891730 

-0.893485 

-0.20 

263.8 

276.0 

-4.6 

7 

2s 2 2p 2 ( 3 P)3s 2 P 

-0.858863 

-0.858604 

0.03 

180.0 

191.2 

-6.2 

8 

2s 2p 4 2 S 

-0.797926 

-0.769659 

3.54 

0.0 

0.0 

0.0 

9 

2s 2 2p 2 ( 3 P)3p 2 S° 

-0.722911 

-0.725836 

-0.40 

0.0 

0.0 

0.0 

10 

2s 2 2p 2 ( 3 P)3p 4 D° 

-0.695903 

-0.698399 

-0.36 

271.7 

287.8 

-5.9 

11 

2s 2 2p 2 ( 1 D)3s 2 D 

-0.695301 

-0.690584 

0.68 

1.0 

-0.2 

121.1 

12 

2s 2 2p 2 ( 3 P)3p 4 P° 

-0.681940 

-0.684113 

-0.32 

138.1 

148.7 

-7.7 

13 

2s 2 2p 2 ( 3 P)3p 2 D° 

-0.652751 

-0.652963 

-0.03 

190.6 

196.3 

-3.0 

14 

2s 2 2p 2 ( 3 P)3p 4 S° 

-0.648007 

-0.649264 

-0.19 

0.0 

0.0 

0.0 

15 

2s 2p 4 2 P 

-0.643687 

-0.621974 

3.37 

168.4 

162.1 

3.8 

16 

2s 2 2p 2 ( 3 P)3p 2 P° 

-0.629299 

-0.628628 

0.11 

59.8 

65.0 

-8.7 

17 

2s 2 2p 2 ( 1 D)3p 2 F° 

-0.496954 

-0.492036 

0.99 

23.6 

23.4 

1.0 

18 

2s 2 2p 2 ( 1 D)3p 2 D° 

-0.485860 

-0.481127 

0.97 

21.4 

15.5 

27.6 

19 

2s 2 2p 2 ( 1 S)3s 2 S 

-0.479897 

-0.466390 

2.81 

0.0 

0.0 

0.0 

20 

2s 2 2p 2 ( 3 P)3d 4 F 

-0.472386 

-0.471490 

0.19 

234.1 

238.9 

-2.0 

21 

2s 2 2p 2 ( 3 P)3d 4 P 

-0.462644 

-0.461378 

0.27 

139.8 

135.0 

3.4 

22 

2s 2 2p 2 ( 1 D)3p 2 P° 

-0.462528 

-0.457352 

1.12 

46.6 

45.8 

1.9 

23 

2s 2 2p 2 ( 3 P)3d 4 D 

-0.460410 

-0.459278 

0.25 

42.2 

38.7 

8.2 

24 

2s 2 2p 2 ( 3 P)3d 2 F 

-0.459102 

-0.458015 

0.24 

162.9 

160.9 

1.3 

25 

2s 2 2p 2 ( 3 P)3d 2 P 

-0.453929 

-0.449836 

0.90 

114.1 

112.9 

1.0 

26 

2s 2 2p 2 ( 3 P)3d 2 D 

-0.445033 

-0.443454 

0.35 

51.8 

48.6 

6.2 


Table 2. The mean energies and total fine-structure splittings of the levels of the configurations ls 2 2s 2 2p 2 ( 3 P j)4f, 5f relative to the 
first ionisation limit. TSTE and TSTT are the experimental and theoretical total fine-structure splitting of each group of levels while 
%S is the percentage difference. Other aspects of the table are as explained in Table [l] 


Configuration 

f'ex 

Eth 

%E 

TSTE 

TSTT 

%S 

2s 2 2p 2 ( 3 P 0 )4f 

-0.251329 

-0.251097 

0.09 

1.4 

1.9 

36.0 

2s 2 2p 2 ( 3 Pi)4f 

-0.250321 

-0.250093 

0.09 

73.5 

70.1 

-4.6 

2s 2 2p 2 ( 3 P 2 )4f 

-0.247971 

-0.247825 

0.06 

230.9 

210.5 

-8.8 

2s 2 2p 2 ( 3 P 0 )5f 

-0.160683 

-0.160573 

0.07 

1.5 

1.9 

27.0 

2s 2 2p 2 ( 3 P!)5f 

-0.159670 

-0.159557 

0.06 

49.1 

46.6 

-5.1 

2s 2 2p 2 ( 3 P 2 )5f 

-0.157735 

-0.157653 

0.05 

111.8 

100.0 

-10.6 
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Table 3. Comparison of level energies in cm -1 within the ls 2 2s 2 2p 2 ( 3 P)4f and ls 2 2s 2 2p 2 ( 3 P)3d configurations between the experimental 
values (E ex ), the work of Liu et al (1995]) (LSBC), and the current work (CW), where the given energies are relative to the lowest level 
in each configuration. The first three columns are reproduced from Table (3) of |Liu et a/||l995| ). 


Level 

-E-ex 

LSBC 

CW 

4f D[3]°5 / 2 

0.0 

0.0 

0.0 

4f D[3]° 7 /2 

1.4 

1.0 

1.9 

4f G[3]° 5/ 2 

66.2 

66.3 

67.4 

4f G[3]° 7 /2 

69.4 

69.5 

71.0 

4f D[2]° 3/ 2 

122.7 

122.7 

124.2 

4f D[2]° 5/2 

123.5 

123.2 

125.4 

4f G[4]°g /2 

137.7 

137.7 

135.2 

4f G[4]° 7/2 

139.6 

139.8 

137.5 

4f D[l]° 3/ 2 

222.3 

221.9 

227.3 

4f D[l]° 1/2 

222.4 

222.4 

226.9 

4f G[5]° 11/2 

287.5 

287.5 

286.8 

4f G[5]°g /2 

293.6 

294.2 

293.7 

4f F[2]° 3/2 

393.7 

393.3 

378.0 

4f F[2]° 5/ 2 

397.8 

397.3 

382.8 

4f F[3]° 7/2 

433.3 

432.9 

417.5 

4f F[3]° 5/2 

435.8 

435.4 

417.5 

4f F[4]° 9/ 2 

446.1 

445.8 

429.5 

4f F[4]° 7/2 

453.3 

452.8 

437.7 

3d 4F 3/ 2 

0.0 

0.0 

0.0 

3d 4 F 

53.9 

54.5 

55.3 

5/2 

3d F 7/ 2 

131.8 

133.0 

134.8 

3d 4 F 

234.1 

236.1 

238.9 

3d 4 P , 

1166.6 

1158.6 

1211.5 

3d 4 P 

1239.8 

1234.2 

1282.3 

3/2 

3d 4p 1/2 

1306.4 

1303.5 

1346.4 

3d 4 D 

1415.5 

1415.2 

1447.6 

. 1/2 

M 4D 3/2 

1449.9 

1449.6 

1481.0 

3d 4 D 

1451.4 

1453.0 

1478.7 

3d 4D 7/ 2 

1457.7 

1460.9 

1486.4 

3d 2 F 

1500.2 

1502.0 

1525.3 

3d 2 F , 

1663.1 

1665.6 

1686.2 

3d 2 P ; 

2134.4 

2131.4 

2488.5 

3d 2 P 

2248.4 

2260.5 

2601.4 

1/2 

3d 2 D , 

3106.7 

3105.4 

3186.6 

3d 2 D 

3158.5 

3161.7 

3235.2 

5/2 
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Table 4. Branching ratios for the strongest 3d-3p transitions in Case B as obtained by |Liu et al\ |l995[ ) (LSBC) and the current work 
(CW) where A is the air wavelength. 


Transition 

A[A] 

LSBC 

CW 

4 F 9 /2 

- 4 d° 

4075.86 

1.000 

1.000 

4F 7/2 

- 4 d°, 

4092.93 

0.119 

0.119 

4 FV/2 

- 4 d° ; 

4072.15 

0.864 

0.862 

4 F 5 /2 

- 4 d° 

4106.02 

0.007 

0.007 

4 F 5 /2 

- 4 d° / 

4085.12 

0.216 

0.216 

4 F 5 /2 

- 4 d° 

4069.89 

0.761 

0.759 

4 F 3 /2 

- 4 d° 

5/2 

4094.14 

0.016 

0.017 

4 F 3 /2 

- 4 d° 

4078.84 

0.264 

0.264 

4 F 3 /2 

- 4 d° 

1/2 

4069.62 

0.715 

0.715 

4 E>7/2 

- 4 d° 

3882.19 

0.137 

0.134 

4 F>7/2 

- 4 d° 

5/2 

3863.50 

0.016 

0.016 

4 E>7/2 

- 2 d° 

4751.28 

0.021 

0.021 

4 F>7/2 

4p° 

5/2 

4119.22 

0.378 

0.370 

4 P 5 /2 

- 4 d°, 

3907.46 

0.054 

0.059 

4 P 5 /2 

4po 

4169.23 

0.153 

0.133 

4 P 5 /2 

- 4 d° ; 

3893.52 

0.013 

0.014 

4 P 5 /2 

4po 

4153.30 

0.449 

0.458 

4 P 5 /2 

4oo 

^3/2 

4924.53 

0.289 

0.283 
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Table 5. Branching ratios for the strongest 4f-3d transitions. All other aspects are as for Table [4] 


Transition 

a[A] 

LSBC 

cw 

G[5]? 1/2 

- 

4 F 9 /2 

4089.29 

1.000 

0.999 

G[5]§ /2 

- 

4 F 9 /2 

4088.27 

0.023 

0.023 

G[5]§ /a 

- 

4 f 7 / 2 

4071.24 

0.204 

0.204 

G[5]§ /a 

- 

2 f 7 / 2 

4342.01 

0.702 

0.695 

G[5]fl /a 

- 

4 D 7 / 2 

4303.61 

0.072 

0.076 


- 

4 F 9 /2 

4114.51 

0.019 

0.017 

GWS /a 

- 

4 f 7 / 2 

4097.26 

0.735 

0.731 

G[4]g /a 

- 

2 f 7/2 

4371.62 

0.125 

0.118 

G[4]g /a 

- 

4 D 7 / 2 

4332.70 

0.121 

0.133 

G[4]? /a 

- 

4 F 7 / 2 

4096.94 

0.046 

0.044 

G[4]? /a 

- 

4 Fs/2 

4083.90 

0.428 

0.423 

G[4]? /a 

- 

2 F 5 /2 

4340.33 

0.354 

0.320 

G[4]? /a 

- 

4 D 5 /2 

4331.17 

0.123 

0.159 

G[4]° 7/2 

- 

2 E> 5/2 

4677.07 

0.046 

0.052 

G[3]? /a 

- 

4 f 7 / 2 

4108.76 

0.054 

0.054 

G[3]? /a 

- 

4 D 7 /2 

4345.56 

0.050 

0.052 

G[3]? /a 

- 

4 F 5 /2 

4095.64 

0.296 

0.280 

G[3]? /a 

- 

2 F 5 /2 

4353.59 

0.162 

0.117 

G[3]? /a 

- 

4 F> 5/2 

4344.38 

0.189 

0.221 

G[3]? /a 

- 

4 P 5 /2 

4291.26 

0.246 

0.273 

4 [4] 9/2 

- 

j f 9 / 2 

4062.93 

0.150 

0.152 

F [4]g/ 2 

- 

4 f 7 / 2 

4046.12 

0.015 

0.019 

F l% 2 

- 

2 f 7 / 2 

4313.44 

0.153 

0.166 

F l% 2 

- 

4 D 7 /2 

4275.55 

0.681 

0.662 

F[4]?/ 2 

- 

4 f 7/2 

4044.95 

0.024 

0.022 

F[4]? /a 

- 

2 f 7/2 

4312.11 

0.098 

0.103 

F[4]? /2 

- 

4 F^ 7 / 2 

4274.25 

0.035 

0.030 

F[4]° /2 

- 

4 F 5 / 2 

4032.24 

0.017 

0.019 

F[4]? /a 

- 

4 F ) 5/2 

4273.10 

0.081 

0.070 

F[4]° /2 

- 

2 D 5 /2 

4609.43 

0.725 

0.733 

F[4]? /2 

- 

4 P 5 /2 

4221.70 

0.017 

0.020 

F[3]° /2 

- 

4 F 9 /2 

4065.05 

0.016 

0.016 

F[3]° /2 

- 

4 F 7 /2 

4048.22 

0.094 

0.099 

F[3]° /2 

- 

2 F 7/2 

4315.83 

0.046 

0.044 

F[3] 7 / 2 

- 

4 D 7 /2 

4277.90 

0.154 

0.159 

F[3] 7/2 

- 

2 F 5 /2 

4285.68 

0.298 

0.354 

F[3] 7/2 

- 

4 F ) 5/2 

4276.75 

0.314 

0.265 

F[3] 7/2 

- 

2 D 5 /2 

4613.68 

0.074 

0.059 

D[3]° /a 

- 

4 D 7 /2 

4358.44 

0.026 

0.021 

D[3]? /2 

- 

4 F 5 /2 

4107.09 

0.173 

0.190 

D[3] 7/2 

- 

2 F 5 /2 

4366.53 

0.050 

0.047 

D[3] 7/2 

- 

4 F ) 5/2 

4357.25 

0.079 

0.091 

D[3] 7/2 

- 

2 D 5 /2 

4707.50 

0.014 

0.016 

D[3] 7/2 

- 

4p 5/2 

4303.82 

0.652 

0.630 
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